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Abstract 

Photonic nanojets have been brought into attention ten years ago for potential application in ul¬ 
tramicroscopy, because of its sub-wavelength resolution that can enhance detection and interaction 
with matter. For these novel applications under development, the optical trapping of a sphere acts 
as an ideal framework to employ photonic nanojets. In the present study, we generated nanojets by 
using a highly focused incident beam, in contrast to traditional plane waves. The method inherits 
the advantage of optical trapping, especially for intracellular applications, with the microsphere 
in equilibrium on the beam propagation axis and positioned arbitrarily in space. Moreover, owing 
to optical scattering forces, when the sphere is in equilibrium, its center shifts with respect to the 
focal point of the incident beam. However, when the system is in stable equilibrium with a con- 
hguration involving optical tweezers, photonic nanojets cannot be formed. To overcome this issue, 
we employed double optical tweezers in an unorthodox configuration involving two collinear and 
co-propagating beams, the precise positioning of which would turn on/off the photonic nanojets, 
thereby improving the applicability of photonic nanojets. 

Keywords: Photonic nanojets. Generalized Lorenz-Mie theory. Beam shape coefficients, Mie 
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1. Introduction 

A photonic nanojet (PNJ) is a narrow (subwavelength) and elongated region with high in¬ 
tensity located at the shadow-side surface of an illuminated loss-less dielectric microcylinder or 
microsphere mm- In the present study, we focus on PNJs generated by a dielectric microsphere 
acting as a focusing lens, which results in the build up of spherically aberrated rays in the fo¬ 
cal region. This build up of rays in the spatially localized high-energy-density region constitute 
the external caustic mm- Therefore, PNJs are not a new phenomenon, but they have attracted 
renewed interest because of recent technological advances allowing the exploitation of this high- 
energy-density region. Nonetheless, contrary to a possible understanding via geometrical optics or 
catastrophe theory, we aim to obtain an exact solution of Maxwell’s equations by employing the 
generalized Lorenz-Mie theory (GLMT) [5], as in the study by Devilez et al. [6]. Such focusing of 
light using a microsphere has also been investigated by Kofler et al. [7] using the uniform caustic 
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asymptotic method, to obtain analytical expressions for the intensity by matching the geometrical- 
optic solutions with Bessoid integrals. However, for optical trapping applications, the scattering 
object is a microsphere of size comparable to the illuminating wavelength and consequently outside 
the regime of geometrical optics. 

Optical trapping results from the change in linear momentum of a beam scattered by a micro¬ 
sphere, producing a resultant restoring force. The most common setup for laser trapping is of opti¬ 
cal tweezers BM- PNJs and optical tweezers (OT) have attracted attention because, in combina¬ 
tion, they can be highly useful in applications including nanoscale processing [iniciiiiiiiiaiiiiiisi, 
high-resolution microscopy [HlIIllIITlIIHlIlg], and enhanced inelastic spectroscopy, such as Raman 
scattering l20l[2ll[22l[23], coherent anti-Stokes Raman scattering |24l [25] and fluorescence [261 [27] , 
or elastic enhancement through backscattering from nanoparticles [28| . A recent study that aimed 
to position PNJs in a controlled manner involved a microsphere attached to a movable micropipette 
[29] and magnetically controlled Janus spheres |3^. In addition, an optically trapped non-spherical 
dielectric particle in a Bessel beam was employed to generate PNJs from another pulsed laser for 
direct laser writing m- The disadvantage of this approach is in the engineering of particular non- 
spherical particles, each exhibiting a characteristic PNJ. Other studies on optical forces and PNJs 
are related to the radiation pressure effect produced by a PNJ on nearby nanoparticles [571 [551 [51] . 
In contrast, the scope of this study is different; we investigate the PNJs generated from an optically 
trapped microsphere in terms of the incoming beam, trading the complexity of structuring special 
shaped particles for the ease of manipulating the trapping beam. 

It is widely known that the main features of PNJs are waists smaller than the diffraction limit 
and propagation over several wavelengths without significant diffraction. These features lead to 
the potential application of PNJs for developing spectroscopic methods with high spatial resolution 
and high detection sensitivity through backscattering enhancement. Therefore, the characteristics 
of a nanojet have been described with a few parameters: the location of the intensity maximum, 
the distance from this maximum to the sphere surface (radial shift), the distance from the intensity 
peak to the point where the intensity decays to times the initial intensity in the direction of 
the beam (decay length), and the PNJ width, which is waist of the jet [35] . The behavior of 
PNJs for a given wavelength is as follows: the PNJ waist widening is proportional to the relative 
refractive index, and the PNJ lengthening is due to the increasing sphere size. This transverse 
confinement (PNJ waist) has been approximated as a Gaussian, while a part of the longitudinal 
confinement (decay length) is approximated as a Lorentzian profile [56|. In the present work, it 
is shown that a bessoidal-type surface best describes the shape of light confined by these PNJs 
originating from the interference between the scattered field and the incident beam |6] and therefore 
the extinction component of the Poynting vector. 

The modeling of PNJs originating from plane-wave incidence involves only four parameters: 
the refractive indices of the particle and its surroundings, wavelength of the incident wave, and 
radius of the sphere. However, with an incident focused beam, which is required for an optical 
trap, the complexity of the modeling increases with additional parameters: the numerical aperture 
and location of the waist with respect to the scatterer. The positioning of the incident focused 
beam with respect to the scatterer is very important, and we demonstrate that it is responsible for 
switching the PNJs on and off for an optically trapped microsphere. 
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2. Theoretical model of photonic nanojets 

Any Maxwellian beam can be expressed by its electromagnetic field in terms of partial waves. 
We start with the general expression for the electromagnetic fields of interest: the incident and 
scattered fields. Later, we simplify this to a two-dimensional problem, which is useful for a 2D plot 
and subsequently for a ID profile such as the longitudinal intensity of the PNJ. 

2.1. Generalized Lorenz-Mie theory 

The basic idea of GLMT is that a beam, which is expressed by a solution to Maxwell’s equations 
(i.e. Maxwellian beam), can be written as an infinite series of spherical functions and spherical 
harmonics, each multiplied by a coefficient called a beam shape coefficient (BSC). These BSCs 
completely describe the incident, internal, and scattered electromagnetic fields in terms of partial- 
waves series. The correct determination of these radially independent BSCs allows for the precise 
determination of the observed electromagnetic phenomena [5j. 

The present notation follows that of previous works [3711381139], expressing the electromagnetic 
fields in terms of spherical vector wave functions even though the original formulation of CLMT 
was in the framework of the Bromwich scalar functions llQlil]. The microsphere scatterer has a 
radius a and real refractive index n, and the incident beam is directed towards the positive 2 ; axis of 
the rectangular coordinate system. The incident wave has a time dependence exp{—iujt) (omitted 
for clarity), wavelength A, and wave number k = 27r/A. The scatterer and the surrounding medium 
are homogeneous, isotropic, and non-magnetic. The incident and scattered electromagnetic fields 
can be written as follows: 


00 p 


= ^0 E E + G™Mp,(r), 

p=l q=-p 

( 1 ) 

00 P 

ll^nc = ^0 E E GprMpg(r) - G™Np,(r), 

p=lq=-p 

( 2 ) 

00 P 

Esca — -^0 ^ ^ ^ ^ Q'pgNpg(r) + 6pg]V[pgf(r), 

P=1 q=-p 

( 3 ) 

00 p 

( 4 ) 


p=lq=-p 


where the terms involving the spherical functions and vectorial spherical harmonics can be 
abbreviated as 


Np(j(r) ~ ^ (^) 

Mpg(r) = Zpq{kr)Xpq{e, 0), (6) 

where Zpq{kr) denotes spherical Bessel or spherical Hankel functions, depending on whether 
it expresses the incident or scattered field, respectively. The vectorial spherical harmonics are 
defined as Xpg(r) = 'LYim{r )/where L)m(r) denotes the scalar spherical harmonics and 
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L = —ir X d/dr is the angular momentum operator in direct space. The scattering partial-wave 
coefficients apq and bpq are related to the Mie scattering coefficients by 


^pq 

m'i/p{mx)il)'p{x) — 'ijjp{mx)ipp{x) 

( 7 ) 

r<TM 

PQ 

Tp'p{mx)Cp{x) - rmpp{mx)Cp{x) ’ 

bpq 

m'ip/{mx)i/p{x) — 'i/ip{mx)il)'p{x) 

( 8 ) 

nTE 

pq 

il}p{mx)i'p{x) - m'i/'p{mx)ip{x) ' 


The minus sign was included to adopt the more frequent notation for the Mie coefficients Op 
and bp. The coefficients Gp^ and Gp^ represent the expansion coefficients (BSCs) of the incident 
fields. The BSCs are defined as follows: 

^TM ^ [\esm6 r\c/Yp\{6,<p)B-r, ( 9 ) 

Bo Jo Jo 

£ dOsmB l\c/Yp*^{9,m-r, (10) 

where gp = kr/{jp{kr)y^n(nj^l)) was introduced to shorten the equation. Note that the 
explicit cancellation of the radial dependence in gp has been the basis of various approximation 
techniques for the BSC um- Recently, the analytical solution for any type of Maxwellian beam 
has been demonstrated [32] . 

For all the incident beams described henceforth, the axis of symmetry coincides with that of the 
sphere, and the scenario is termed on-axis. This is due to the fact that the equilibrium position for 
a homogeneous sphere lies on the axis of symmetry. For this particular conhguration, the double 
sum in the partial wave expansion, Eqs. ([T])-©, simplifies because of contributions exclusively 
from g = ±1 [431 [39]. For this case, it would be best to introduce a new vectorial function to 
describe our electromagnetic fields: 

N±(r) = Np,i(r)±Np,_i(r), (11) 


M±(r) = Mp,i(r)±Mp,_i(r), 


( 12 ) 


thereby simplifying the analytical solution up to this point using the symmetrical relations 
for q = ±1. Further simplification of the BSC would require knowledge of the exact vectorial 


electromagnetic fields near the focus, as will be shown in section 2.4 


2.2. Plane wave incidence 

Any incident field can be expressed as an infinite series of vectorial spherical harmonics. As an 
initial approach, an arbitrary polarized plane wave travelling along the z-axis (k = kz) is used to 
validate the results and compare them with known results. The electric field is 


B = EoP^-^ {Px,Py,0), 


(13) 


where the polarization is described by the components of the Jones vector, {px,Py)- The 
magnetic field can be easily determined from the relation 
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(14) 


The radial components of the fields in spherical coordinates are 


H = — X E = i?oe* 


E • f = Eqb 


ikr cos 9 


H • r = Hoe 


ikr cos 9 


sin 6 {px cos (j) + Py sin <p), 
sin 6 {px sin 4> ~ Py cos (f )), 


(15) 

(16) 


where we used f = sin0cos(/)x + sin0sin(/)y + cosOz. The integrals over the solid angle are 
presented in Appendix A, and it can be seen that only the q = ±1 components are present. The 
plane-wave (pw) BSCs are 


G 


{^ipx - Py ), ip^ T iPy ), 


TE,pw _ 


(17) 


GfT = v'vr(2p+l)i^’. (18) 

2.3. Description of highly focused beam 

The commonly employed Davis description of Gaussian beams, a perturbation method that 
introduces higher-order corrections to the paraxial approximation, yields a condition that is rarely 
met for strong focusing (as in optical tweezers). Consequently there is a need to go beyond a 
paraxial-expansion description. To achieve this aim, we resort to an early description based on the 
angular spectrum representation, which is later applied to optical trapping to result in a solution 
for a completely arbitrary vectorial Maxwellian beam 1371. This new approach has the benefit of 
providing an analytical expression for the BSCs as a function of beam position and polarization 
directly with respect to the sphere without the need to rely on vector translation theorems requiring 
high computational cost. Results obtained by using this approach have recently been confirmed 
independently through an aberrated photothermal microscopy measurement to be in excellent 
agreement with the theory |44j . 

To present a complete study of PNJs in optical tweezers, we start with a description of the 
highly focused incident beam in the framework of GLMT. The on-axis case is considered here, since 
the equilibrium trap position for a homogeneous sphere lies on the propagation axis. The exact 
BSCs for an on-axis highly focused beam (fb) positioned at zq with respect to the sphere center 
are [39] 


r'pw 

GfJ’ = ikfE^G Gp 


f 


[Tr^ia) + T^{a)] 


P{P + 1) 

In terms of the two angular functions, which involves associated Legendre functions, 


( 20 ) 


4(0) = 


Pp{cos9) 

sin0 


r’(f) = 


( 21 ) 


5 







2.4- Approximations for fields along the incident plane (2D) 

When plotting the fields in 2D, in the x — z plane, we will take (j) = 0 but 0 7 ^ 0 in Eqs. 
(|ll[)-p!2|). The vectorial incident and scattered fields from Eqs. 0 -@ becomes for the on-axis 
case, 
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For the plane-wave BSCs, (i.e. Gp = G^), results in 
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Note that for the Poynting vector in the x — z plane, we only need the r and 9 components. 
We are thus interested in the following products from the vectorial fields above: 


Srf + See = f {EeH^ - E^Hg) + 9 {E^Hr - ErH^). 


(30) 


For the particular 2D case, since the photonic nanojet is extremely close to the z-axis immedi¬ 
ately outside the spherical scatterer, the small-angle approximation for small 9 is justified: 


f’.ncosP)» (-i).^/|±|y+:j,(,+(^p). 

plicose) = -\/p(p + 1) V Ji( Vp(p + 1)0). 

V sm (7 


(31) 

(32) 


This small-angle approximation yields numerically faster 2D plots that are not noticeably dif¬ 
ferent from the plot obtained through the use of the exact expression of the associated Legendre 
functions. 


2.5. Approximations for fields along the z-axis (ID) 

When plotting only the values of the longitudinal intensity profile of the PNJ, we are interested 
only in the fields along the z-axis (i.e., 0 = 0); the angular functions (Eqs. ([^-Q) for 0 —0 are 
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de 
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0-5-0 2 


This greatly simplifies the field, especially for a fast calculation of the longitudinal profile using 
the Poynting vector of section 2.6, resulting in fields with no r-component. 


2.6. Time-averaged Poynting veetor 

By examining the time-averaged Poynting vector, we aim to elucidate the vectorial character¬ 
istics of the PNJ and its dependence on the incident polarization. Previous results [l5] considered 
the total field for scattering by a microsphere to represent the PNJ. Here, emphasis is given to 
the extinction component of the Poynting vector [46] as the origin of the PNJ, as suggested in [6|. 
We can now determine the Poynting vector for the incident, scattered, and interference fields. The 
electromagnetic fields outside the spherical scatterer are 


Etot — Tjinc T Escaj T^tot — T^inc T H^ca- 
The time-averaged Poynting vectors are as follows: 

{S^nc) = l^(E^ncXU*^,), 

(S.,a) = 25?(E.caXH:,J, 

(Sexi) = (Etot X Ul,) = (E,„e X H:,J + (E,,, X H* 


(37) 

(38) 

(39) 

(40) 


where Sine, Ssca, and Sext are, respectively, the incident, scattered, and interference (between the 
scattered and incident beams) Poynting vectors. We now carefully examine each Poynting vector 
by using Gp = -y7r(2/r+T)iP for the case of plane waves. Since by definition, \px\^ + \py\‘^ = 1, we 
recover an on-axis intensity profile independent of the chosen incident polarization as follows: 
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• z 






H 2 


P=1 


d{kz) 


(41) 


By using the relations and identities of spherical functions, we recover the Poynting vector for 
the incident field as 


(Sine • i) = 


^e 


ikz 


(42) 


This is expected, as the sum term must converge in Eq. (41) to a particular value because the 
incident plane wave is properly defined as 

{Sine) = {EoHq{pxX +Pyy) x (-p*x -(-p*y)) = {EqHqz) . 
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(43) 












This plane-wave analysis validates the present equations. For Sgca, we start with 


(Ssca • 2 ) = Op +ibph^p \kz) 

^(2p-M)(-i)P b* ^(p + 1) ^^ - h^pllikz)^ - ia*ph':^'>*{kz) (44) 

It is tempting to approximate the spherical Hankel function in its asymptotic form, since the helds 
are in the region outside the spherical scatterer, i.e., kr > ka. For a low p-index, kr » p. 
However, when p increases, the Mie coefficient decreases; in such a case, we expect the asymptotic 
form to be a very good approximation. Unfortunately the asymptotic expansion does not represent 
the Poynting vector up to twice the size factor. 

Finally, we are interested in determining the Poynting vector related to the interference between 


the incident and scattered beams. From Eq. (40), this would be the real part of a linear combination 


of the spherical Hankel function with a plane wave, which is described in the form of the equation 
below: 


(Sext ■ z) = -?R. 


■^ikz Y^(^-i)Ph(^>{kz)Bp - ie-^^^Y,^Ph^^\kz)Ap 

p=0 p=0 


(45) 


Note that Ap and Bp are linear combinations of the complex Mie coefficients Op and bp, and the 
position dependence of the Poynting vector originates from the terms containing (kz) in Eq. (45). 


2. 7. Optical forces 

Once the BSCs are calculated, the force components on an optically trapped microsphere are 
determined via the Maxwell stress tensor formalism, providing new insights owing to the analytical 
nature of the solution. Einally, to represent the PNJ due to an optically trapped microsphere, we 
first determine the location along the propagation axis where stable equilibrium is verified. The 
position is determined from an optical-force curve with respect to the trap position from the 
microsphere center. This type of force prohle has been applied previously in the investigation of 
morphology dependent resonances of an optically trapped microsphere [38], as well as locations 
along the axis [39|. At this stable position, an energy-density plot can be drawn to identify the 
existence/location of PNJs for different optical trapping parameters (wavelength, refractive indexes, 
polarization, and numerical aperture). 

We know that the optical force along the z-axis is 
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For the on-axis case, only p = ±1 remains: 
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The cross products of BSCs for the focused beam can be simplified to 
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This expression for the axial component of the force must be determined for each position of 
the beam focus with respect to the sphere center, which involves determining the BSCs for each of 
these positions. Owing to the linearity of the electromagnetic fields, we can add another focused 
field simply by adding BSCs. Let us denote, zi and Z 2 as the focal positions of the two independent 


focused beams, instead of the previous notation of zq for the single beam in Eq. (20). For the 
double beam the corresponding BSC is 


/^pw fOmax 


p{p + 1) 


rGX-max / \ o 

/ dasinaVcOS^ ^ ^-{fsinaM ^ ^ 

(52) 

which can be simplified if the beam positions are symmetric (i.e., when zi = zq + d/2 and 
Z 2 = zo — d/2, where d represents the spaceing between both focus) as follows: 


^pw pOimax ^ 

Gl^ = 2ikfF^f J dasina\/bd^e*^^o™"“cos(A:d/2 cosa)e-(^"““/‘^) [irlia) + T^{a)] 

^ ^ ° (53) 


3. Results and discussion 

Mathematica (version 10, Wolfram Inc.) was chosen as the computational software for perform¬ 
ing the numerical simulation presented in this section because of its arbitrary-precision computa¬ 
tion, error tracking, and numerical libraries. A suitable truncation number is chosen to terminate 
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the infinite sum, according to m, and the floating-point machine error is chosen. Note that this is 
correct only for the total cross section, as has been emphasized recently |48j . However, for specific 
fields instead of cross sections, a similar rule is applicable in which the maximum number of terms 
greater than the largest kr in the plot region should be chosen. 


3.1. PNJ as a Bessoidal-type surface 

Starting with conventional plane-wave incidence, we determine the Poynting vector of the in¬ 
terference between the incident and scattered fields in a region immediately after the scatterer for 
a 2-/im (Figureand 4-/rm microsphere (Figure]^. 


'? 
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02 
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Figure 1: Contour and vector plot of the intensity and time-averaged Poynting vector, respectively 
for the interference between the incident and scattered field. The result is for a 2-/xm polystyrene 
microsphere (n=1.59) in water (n=1.33) under plane-wave incidence at A=633 nm. 



Figure 2: Contour and vector plot of the intensity and time-averaged Poynting vector, respectively 
for the interference between the incident and scattered fields. The result is for a 4-/im polystyrene 
microsphere (n=1.59) in water (n=1.33) under plane-wave incidence at A=633 nm. 


The htting function for the transverse and longitudinal intensity profiles of a PNJ has been 
described as a Gaussian and a Lorentzian, respectively [36], noting that the intensity distribution 
of photonic jets along the z-axis (longitudinal direction) is not symmetric with respect to the 
maximum intensity. From the results presented here (Figure]^, we can obtain a dominant function 
that adequately represents the profiles of the PNJ. The transverse profile is better described by 
spherical (from Eq. (|45[)) and cylindrical Bessel functions (from the small-angle approximation, 
Eq.(32)), as can be observed by the htting (Eigurej^. 
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(a) Longitudinal intensity profile (b) Transverse intensity profile 


Figure 3: PNJ intensity profile for plane-wave incidence. 
Figure at the maximum intensity. 


The result is for the 2-/rm sphere of 



Figure 4: Fitting of the PNJ transverse profile using a cylindrical Bessel type versus that using a 
Gaussian function. The result is for the 2-/im sphere of Figure |3bl clearly showing a better fitting 
with the Bessel-type function. 


3.2. PNJ from a focused beam 

A PNJ is only observed when the focus is close to the sphere, as has been pointed out by 
Lecler et al. [l5], and the behaviour of the PNJ with respect to the focal position of the beam was 
investigated by Devilez et al. [6]. Therefore, for a highly focused beam at the sphere center, which 
is close to the stable trapping configuration of an optical tweezer, we obtain no PNJ as illustrated 
in Figure 

To observe any appreciable PNJ, we have to shift the incident-beam focus close to the surface. 
By placing it 2 fj,m from the sphere center, we observed the PNJ of Figure Moreover, we 
examined the PNJ transverse profile as a function of beam position, as shown in Figure It can 
be observed that the PNJ due to the interference of the incident and scattered fields, inherits the 
intensity-decay characteristic of the scattered field and localization from the highest-intensity 
position of the incident field. 

The challenge here is to obtain a PNJ from an optically trapped sphere. In an optical tweezer, 
we should not have a PNJ, since the microsphere center is brought near the focus of the beam, in 
a configuration similar to that of Figure 
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Figure 5: Contour and vector plot of the intensity and time-averaged Poynting vector, respectively 
for the interference between the incident and scattered fields. The result is for the 2-/rm sphere 
of Figure 3b, but with a highly focused [NA=\.2f)) Gaussian beam (waist=2.5mm before at the 


objective back aperture) placed at the origin {zq = O/rm). 



Figure 6: Contour and vector plot of the intensity and time-averaged Poynting vector for the 
interference between the incident and scattered fields. The result is for the 2-|Um sphere of figure 


3b, but with a highly focused {NA=1.2b) Gaussian beam (waist=2.5mm before at the objective 


back aperture) located outside the sphere {zq = 2^m). 


3.3. Optical forces in a microsphere 

The optical force profile of conventional optical tweezers is presented in Figure Note that 
the equilibrium position for the microsphere is located immediately after (~ 0.2 /xm) the highly 
focused beam position, which is due to the scattering forces as the focused beam close to the sphere 
center yields no PNJ. 

To circumvent the absence of a PNJ in conventional optical tweezers, we adopt double optical 
tweezers: not the configuration commonly employed in which each beam axis is parallel to each 
other [39], but in a configuration in which both beams are on the same axis. This leads to the 
assumption of the original counter-propagating traps of 1970 [50|, but here, the two traps are co¬ 
propagating. In this double co-propagating optical-tweezers configuration, an equilibrium position 
will exist if the gradient force from one tweezer exceeds the scattering force from the other. 
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Figure 7: Overlapping plot of the on-axis intensity of the time-averaged Poynting vector for the 
interference between the incident and scattered fields. The result is for the 2-/rm sphere of figure 


3b, but with a highly focused {NA=1.2b) Gaussian beam (waist=2.5mm before at the objective 
back aperture) located at specific positions (zq = Zfocus), he., different radial shifts. The dashed 
line is an fitting though the maximum profile points. 



Figure 8: Axial optical force as a function of microsphere displacement z{^m) from the focus. 
The inset illustrates equilibrium configuration, the beam propagation direction (arrow), the focus 
(plus) and microsphere (dot). 


For the 2-;um sphere and two equal-intensity traps located 2 /rm apart (Figure 9a), because 
of the scattering force from the first trap, the equilibrium position is displaced to about 1.1 /rm 
in the beam-propagation direction. There is a cost associated with using these co-propagating 
traps: the reduction of the trap stiffness. The loss of stiffness for the present case is approximately 
40% with respect to a single optical tweezer. Consequently, the second focused beam is positioned 
within the trapped sphere, and no PNJ is generated (Figure [Ob] ). On increasing the displacement 
between the two traps to 3 //m, we obtain two equilibrium positions, as illustrated in Figure [9^ For 
the microsphere located in the first equilibrium position of z=-1.2 /rm, a focused beam is present 
immediately outside the microsphere, resulting in a PNJ (Figure [95] ). For the second equilibrium 
position, there is no PNJ, as expected. Therefore depending on where the microsphere is located 
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(a) Axial optical-force profile 



(b) Contour intensity plot 



(c) Axial optical force profile Contour intensity plot 


Figure 9: Double co-propagating optical traps: one located at z=-1.0 /rm and z=1.0 ^m, corre¬ 
sponding to (a)-(b), and the other located at z=-1.5 and z=1.5 /rm, shown by (c)-(d). The inset 
illustrates equilibrium configuration, the beam propagation direction (arrow), the focus (plus) and 
microsphere (dot). 


within the two possible trap locations, the PNJ can easily be switched on or off by modulating the 
intensity of both traps. 

4. Conclusion 

In summary, we presented a detailed theoretical investigation of the formal relationship between 
the optical force on a microsphere and its photonic nanojet, which results from the interference 
between the scattered and incident fields. This approach is general, and a complex refractive 
index can be implemented to take into account light absorption by the scatterer. Even though 
double optical tweezers were employed, a similar principle holds as for other beam conhgurations, 
especially those generated by holographic traps. 

Consequently, an understanding between the interaction of an optically trapped microsphere 
and its PNJ field structure is of fundamental importance in optical physics and has practical 
significance in applications such as imaging, nano-lithography, detection, metrology, biophotonics, 
and spectroscopy. Such an understanding also provides an opportunity to manipulate the intensities 
of two traps differently in order to achieve a controlled PNJ near an optically trapped microsphere. 
Moreover, the loss in trap stiffness can be easily overcome using anti-reflection coatings to enhance 
the trapping force m- In the case of trapping soft spheres, such as water drops, as has been 
recently reported [52], laser-induced surface stress would break the spherical symmetry of the 
scatterer leading to spheroidal particles. The PNJs generated by such scatterers have recently 
been reported in Finally when whispering-gallery modes or pulsed beams are of interest, they 
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can also be detected in an optical trap [l9], thereby validating recent simulations on photonic-jet 
shaping [Mj and temporal dynamics of the jet [55]. 

Acknowledgments 

The present work received support from CNPq (308627/2012-1), Conselho Nacional de Desen- 
volvimento Cientifico e Tecnologico and FAPESP (2014/07191-9), Fundagao de Amparo a Pesquisa 
do Estado de Sao Paulo, Brazil. 


Appendix A. Shape-coefficient integrals of the beam 


In this appendix, the analytical integrals for the integration of the BSCs are presented. The 
first integral, in terms of the azimuthal angle (j) is of the type 
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which does not depend on 9. Therefore the remaining integral, in terms of the polar angle 6 
for the plane-wave case, is of the type 


d0sin2 0p9(cos0)e*^^“"'^. (A.2) 

This integral resembles a formerly reported one |56j . the solution of which is 

(i6<sin0P^(cos0)e*^^™®"“®® Jg(/cr sinasin0) = 2i^~^ P^{cos a)jp{kr). (A.3) 

Now, we are interested in taking the limit a —)• 0 of the above integral. Because of the properties 
of the Bessel function and associated Legendre functions, we obtain, for q = +1 and g = — 1, 




r dOsin^ 0Pp^(cos0)e*^"“"® = 2p{p + , 

Iq kr 
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JO - kr 

thereby eliminating the radial-dependency term, jp{kr)/kr, from the BSCs. 


(A.4) 

(A.5) 
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